Harnessing the druggability at orthosteric and allosteric sites of PD-1 for small molecule discovery by an integrated in silico pipeline.

The PD-1/PD-L1 interaction is a promising target for small molecule inhibitors in cancer immunotherapy, but targeting this interface has been challenging. While efforts have been made to identify compounds that target the orthosteric sites, no reports have explored the potential of small molecules to target the allosteric region of PD-1. Therefore, our study aims to establish a pipeline to identify small molecules that can effectively bind to either the orthosteric or allosteric pockets of PD-1. We categorized the PD-1 interface into two hot-spot zones (P-and N-zones) based on extensive analysis of its structural, dynamical, and energetic properties. These zones correspond to the orthosteric and allosteric PPI sites, respectively, targeted by monoclonal antibodies. We used a guided virtual screening workflow to identify hits from ∼7 million compounds library, which were then clustered based on structural similarity and assessed by interaction fingerprinting. The selective and diverse chemical representatives were subjected to MD simulations and binding energetics calculations to filter out false positives and identify actual binders. Binding poses metadynamics calculations confirmed the stability of the final hits in the pocket. This study emphasizes the need for an integrated pipeline that uses molecular dynamics simulations and binding energetics to identify potential binders for the dynamic PD-1/PD-L1 interface, due to the lack of small molecule co-crystals. Only a few potential binders were discovered from a large pool of molecules targeting both the allosteric and orthosteric zones. Our results suggest that the allosteric site has more potential than the orthosteric site for inhibitor design. The identified "computational hits" hold potential as starting points for in vitro evaluations followed by hit-to-lead optimization. Overall, this study represents an effort to establish a computational pipeline for exploring and enriching both the allosteric and orthosteric sites of PPI interfaces, "a tough but indispensable nut to crack".

Non-oncology drug (meticrane) shows anti-cancer ability in synergy with epigenetic inhibitors and appears to be involved passively in targeting cancer cells.

Emerging evidence suggests that chemotherapeutic agents and targeted anticancer drugs have serious side effects on the healthy cells/tissues of the patient. To overcome this, the use of non-oncology drugs as potential cancer therapies has been gaining momentum. Herein, we investigated one non-oncology drug named meticrane (a thiazide diuretic used to treat essential hypertension), which has been reported to indescribably improve the therapeutic efficacy of anti-CTLA4 in mice with AB1 HA tumors. In our hypothesis-driven study, we tested anti-cancer potential meticrane in hematological malignance (leukemia and multiple myeloma) and liver cancer cell lines. Our analysis showed that: 1) Meticrane induced alteration in the cell viability and proliferation in leukemia cells (Jurkat and K562 cells) and liver cancer (SK-hep-1), however, no evidence of apoptosis was detectable. 2) Meticrane showed additive/synergistic effects with epigenetic inhibitors (DNMT1/5AC, HDACs/CUDC-101 and HDAC6/ACY1215). 3) A genome-wide transcriptional analysis showed that meticrane treatment induces changes in the expression of genes associated with non-cancer associated pathways. Of importance, differentially expressed genes showed favorable correlation with the survival-related genes in the cancer genome. 4) We also performed molecular docking analysis and found considerable binding affinity scores of meticrane against PD-L1, TIM-3, CD73, and HDACs. Additionally, we tested its suitability for immunotherapy against cancers, but meticrane showed no response to the cytotoxicity of cytokine-induced killer (CIK) cells. To our knowledge, our study is the first attempt to identify and experimentally confirm the anti-cancer potential of meticrane, being also the first to test the suitability of any non-oncology drug in CIK cell therapy. Beyond that, we have expressed some concerns confronted during testing meticrane that also apply to other non-oncology drugs when considered for future clinical or preclinical purposes. Taken together, meticrane is involved in some anticancer pathways that are passively targeting cancer cells and may be considered as compatible with epigenetic inhibitors.

Template Entrance Channel as Possible Allosteric Inhibition and Resistance Site for Quinolines Tricyclic Derivatives in RNA Dependent RNA Polymerase of Bovine Viral Diarrhea Virus.

The development of potent non-nucleoside inhibitors (NNIs) could be an alternate strategy to combating infectious bovine viral diarrhea virus (BVDV), other than the traditional vaccination. RNA-dependent RNA polymerase (RdRp) is an essential enzyme for viral replication; therefore, it is one of the primary targets for countermeasures against infectious diseases. The reported NNIs, belonging to the classes of quinolines (2h: imidazo[4,5-g]quinolines and 5m: pyrido[2,3-g] quinoxalines), displayed activity in cell-based and enzyme-based assays. Nevertheless, the RdRp binding site and microscopic mechanistic action are still elusive, and can be explored at a molecular level. Here, we employed a varied computational arsenal, including conventional and accelerated methods, to identify quinoline compounds' most likely binding sites. Our study revealed A392 and I261 as the mutations that can render RdRp resistant against quinoline compounds. In particular, for ligand 2h, mutation of A392E is the most probable mutation. The loop L1 and linker of the fingertip is recognized as a pivotal structural determinant for the stability and escape of quinoline compounds. Overall, this work demonstrates that the quinoline inhibitors bind at the template entrance channel, which is governed by conformational dynamics of interactions with loops and linker residues, and reveals structural and mechanistic insights into inhibition phenomena, for the discovery of improved antivirals.

Deciphering the Structural Determinants Critical in Attaining the FXR Partial Agonism.

Elucidation of structural determinants is pivotal for structure-based drug discovery. The Farnesoid X receptor (FXR) is a proven target for NASH; however, its full agonism causes certain clinical complications. Therefore, partial agonism (PA) appears as a viable alternative for improved therapeutics. Since the agonist and PA both share the same binding site, i.e., ligand-binding pocket (LBP), which is highly dynamic and has synergy with the substrate binding site, the selective designing of PA is challenging. The identification of structural and conformational determinants is critical for PA compared with an agonist. Furthermore, the mechanism by which PA modulates the structural dynamics of FXR at the residue level, a prerequisite for PA designing, is still elusive. Here, by using ∼4.5 µs of MD simulations and residue-wise communication network analysis, we identified the structural regions which are flexible with PA but frozen with an agonist. Also, the network analysis identified the considerable changes between an agonist and PA in biologically essential zones of FXR such as helix H10/H11 and loop L:H11/H12, which lead to the modulation of synergy between LBP and the substrate binding site. Furthermore, the thermodynamic profiling suggested the methionine residues, mainly M328, M365, and M450, seem to be responsible for the recruitment of PA. The other residues I357, Y361, L465, F308, Q316, and K321 are also identified, exclusively interacting with PA. This study offers novel structural and mechanistic insights that are critical for FXR targeted drug discovery for PA designing.

Discovery of non-nucleoside oxindole derivatives as potent inhibitors against dengue RNA-dependent RNA polymerase.

A series of thiazole linked Oxindole-5-Sulfonamide (OSA) derivatives were designed as inhibitors of RNA-dependent RNA polymerase (RdRp) activity of Dengue virus. These were synthesized and then evaluated for their efficacy in ex-vivo virus replication assay using human cell lines. Among 20 primary compounds in the series, OSA-15 was identified as a hit. A series of analogues were synthesized by replacing the difluoro benzyl group of OSA-15 with different substituted benzyl groups. The efficacy of OSA-15derivatives was less than that of the parent compound, except OSA-15-17, which has shown improved efficacy than OSA-15. The further optimization was carried out by adding dimethyl (DM) groups to both the sulfonamide and oxindole NH's to produce OSA-15-DM and OSA-15-17-DM. These two compounds were showing no detectable cytotoxicity and the latter was more efficacious. Further, both these compounds were tested for inhibition in all the serotypes of the Dengue virus using an ex-vivo assay. The EC50 of OSA-15-17-DM was observed in a low micromolar range between 2.5 and 5.0 µg/ml. Computation docking and molecular dynamics simulation studies confirmed the binding of identified hits to DENV RdRp. OSA15-17-DM blocks the RNA entrance and elongation site for their biological activity with high binding affinity. Overall, the identified oxindole derivatives are novel compounds that can inhibit Dengue replication, working as non-nucleoside inhibitors (NNI) to explore as anti-viral RdRp activity.

Traversing through the Dynamic Protein-Protein Interaction Landscape and Conformational Plasticity of PD-1 for Small-Molecule Discovery.

Monoclonal antibodies (mAbs) blocking the PD-1/PD-L1 interface have shown remarkable success in treating malignancies, but they may also initiate lethal immune-related adverse events. Small molecules may circumvent the mAb limitations; however, none has entered clinical trials targeting PD-1. Its complex protein-protein interaction interfaces necessitate an atomic-level understanding of recognition and binding mechanisms. Hence, we have aimed to highlight the PD-1's sequence-structure-dynamic-function link with its cognate ligands and diversely reported inhibitors. We focus primarily on the anti-PD-1 mAbs, their mode of actions, and interactions with PD-1 epitopes. The comparison of co-crystals showed that these ligands/inhibitors harness the PD-1's conformational plasticity and structural determinants differentially. The relationship between modulator binding patterns and biological activity is demonstrated using interaction fingerprinting of all reported human PD-1 co-crystals. The significant dynamical events and hot-spot residues underpinned from crystallographic wealth and computational studies have been highlighted to expedite small-molecule discovery.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations.

The main protease, Mpro/3CLpro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the Mpro mechanism of action. A comparative analysis of both forms of Mpro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of Mpro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of Mpro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent Mpro inhibitors.Communicated by Ramaswamy H. Sarma.

Characterizing (un)binding mechanism of USP7 inhibitors to unravel the cause of enhanced binding potencies at allosteric checkpoint.

The ability to predict the intricate mechanistic behavior of ligands and associated structural determinants during protein-ligand (un)binding is of great practical importance in drug discovery. Ubiquitin specific protease-7 (USP7) is a newly emerging attractive cancer therapeutic target with bound allosteric inhibitors. However, none of the inhibitors have reached clinical trials, allowing opportunities to examine every aspect of allosteric modulation. The crystallographic insights reveal that these inhibitors have common properties such as chemical scaffolds, binding site and interaction fingerprinting. However, they still possess a broader range of binding potencies, ranging from 22 nM to 1,300 nM. Hence, it becomes more critical to decipher the structural determinants guiding the enhanced binding potency of the inhibitors. In this regard, we elucidated the atomic-level insights from both interacting partners, that is, protein-ligand perspective, and established the structure-activity link between USP7 inhibitors by using classical and advanced molecular dynamics simulations combined with linear interaction energy and molecular mechanics-Poisson Boltzmann surface area. We revealed the inhibitor potency differences by examining the contributions of chemical moieties and USP7 residues, the involvement of water-mediated interactions, and the thermodynamic landscape alterations. Additionally, the dissociation profiles aided in the establishment of a correlation between experimental potencies and structural determinants. Our study demonstrates the critical role of blocking loop 1 in allosteric inhibition and enhanced binding affinity. Comprehensively, our findings provide a constructive expansion of experimental outcomes and show the basis for varying binding potency using in-silico approaches. We expect this atomistic approach to be useful for effective drug design.

Targeting cryptic-orthosteric site of PD-L1 for inhibitor identification using structure-guided approach.

Approved mAbs that block the protein-protein interaction (PPI) interface of the PD-1/PD-L1 immune checkpoint axis have led to significant improvements in cancer treatment. Despite having drawbacks of mAbs only few a compounds are reported till date against this axis. Inhibiting PPIs using small molecules has emerged as a significant therapeutic opportunity, demanding for the identification of drug-like molecules at an accelerated pace under the hit-to-lead campaigns. Due to the PD-L1's cross-talk with PD-1/CD80 and its overexpression on cancer cells, as well as the availability of its crystal structures with small molecules, it is an enticing therapeutic target for structure-assisted small molecule design. Furthermore, the selection of chemical databases enriched with focused designing for PPI interfaces is crucial. Therefore, in this study we have utilized the Asinex signature library for structure-assisted virtual screening to find the potential PD-L1 inhibitors by targeting the cryptic PD-L1 interface, followed by induced fit docking for pose refinements in the pocket. The obtained hits were then subjected to interaction fingerprinting and ligand-based drug-likeness investigations in order to evaluate and analyze their drug-like qualities (ADME). Twelve compounds qualified for molecular dynamics simulations, followed by thermodynamic calculations for evaluation of their stability and energetics inside the pocket. Two novel compounds with different chemical moieties have been identified that are consistent throughout the simulation, mimicking the interactions and binding energies with BMS-1166. These compounds appear as potential therapeutic candidates to be explored experimentally, thereby paving the way for the development of novel leads as immunomodulators.

Molecular Dynamics Simulations Reveal the Interaction Fingerprint of Remdesivir Triphosphate Pivotal in Allosteric Regulation of SARS-CoV-2 RdRp.

The COVID-19 pandemic has now strengthened its hold on human health and coronavirus' lethal existence does not seem to be going away soon. In this regard, the optimization of reported information for understanding the mechanistic insights that facilitate the discovery towards new therapeutics is an unmet need. Remdesivir (RDV) is established to inhibit RNA-dependent RNA polymerase (RdRp) in distinct viral families including Ebola and SARS-CoV-2. Therefore, its derivatives have the potential to become a broad-spectrum antiviral agent effective against many other RNA viruses. In this study, we performed comparative analysis of RDV, RMP (RDV monophosphate), and RTP (RDV triphosphate) to undermine the inhibition mechanism caused by RTP as it is a metabolically active form of RDV. The MD results indicated that RTP rearranges itself from its initial RMP-pose at the catalytic site towards NTP entry site, however, RMP stays at the catalytic site. The thermodynamic profiling and free-energy analysis revealed that a stable pose of RTP at NTP entrance site seems critical to modulate the inhibition as its binding strength improved more than its initial RMP-pose obtained from docking at the catalytic site. We found that RTP not only occupies the residues K545, R553, and R555, essential to escorting NTP towards the catalytic site, but also interacts with other residues D618, P620, K621, R624, K798, and R836 that contribute significantly to its stability. From the interaction fingerprinting it is revealed that the RTP interact with basic and conserved residues that are detrimental for the RdRp activity, therefore it possibly perturbed the catalytic site and blocked the NTP entrance site considerably. Overall, we are highlighting the RTP binding pose and key residues that render the SARS-CoV-2 RdRp inactive, paving crucial insights towards the discovery of potent inhibitors.

Conformational Characterization of the Co-Activator Binding Site Revealed the Mechanism to Achieve the Bioactive State of FXR.

FXR bioactive states are responsible for the regulation of metabolic pathways, which are modulated by agonists and co-activators. The synergy between agonist binding and 'co-activator' recruitment is highly conformationally driven. The characterization of conformational dynamics is essential for mechanistic and therapeutic understanding. To shed light on the conformational ensembles, dynamics, and structural determinants that govern the activation process of FXR, molecular dynamic (MD) simulation is employed. Atomic insights into the ligand binding domain (LBD) of FXR revealed significant differences in inter/intra molecular bonding patterns, leading to structural anomalies in different systems of FXR. The sole presence of an agonist or 'co-activator' fails to achieve the essential bioactive conformation of FXR. However, the presence of both establishes the bioactive conformation of FXR as they modulate the internal wiring of key residues that coordinate allosteric structural transitions and their activity. We provide a precise description of critical residue positioning during conformational changes that elucidate the synergy between its binding partners to achieve an FXR activation state. Our study offers insights into the associated modulation occurring in FXR at bound and unbound forms. Thereafter, we also identified hot-spots that are critical to arrest the activation mechanism of FXR that would be helpful for the rational design of its agonists.

Interplay among Structural Stability, Plasticity, and Energetics Determined by Conformational Attuning of Flexible Loops in PD-1.

The dynamics and plasticity of the PD-1/PD-L1 axis are the bottlenecks for the discovery of small-molecule antagonists to perturb this interaction interface significantly. Understanding the process of this protein-protein interaction (PPI) is of fundamental biological interest in structure-based drug designing. Food and Drug Administration (FDA)-approved anti-PD-1 monoclonal antibodies (mAbs) are the first-in-class with distinct binding modes to access this axis clinically; however, their mechanistic aspects remain elusive. Here, we have unveiled the interactive interfaces with PD-L1 and mAbs to investigate the native plasticity of PD-1 at global (structural and dynamical) and local (residue side-chain orientations) levels. We found that the structural stability and coordinated Cα movements are increased in the presence of PD-1's binding partners. The rigorous analysis of these PPIs using computational biophysical approaches revealed PD-1's intrinsic plasticity, its concerted loops' movement (BC, FG, and CC'), distal side-chain motions, and the thermodynamic landscape, which are perturbed remarkably from its unbound to bound states. Based on intra-/inter-residues' contact networks and energetics, the hot-spots have been identified that were found to be essential to arrest the dynamical motions of PD-1 significantly for the rational design of therapeutic agents by mimicking the mAbs mechanism.

Identification of potential molecules against COVID-19 main protease through structure-guided virtual screening approach.

The pandemic caused by novel coronavirus disease 2019 (COVID-19) infecting millions of populations worldwide and counting, has demanded quick and potential therapeutic strategies. Current approved drugs or molecules under clinical trials can be a good pool for repurposing through in-silico techniques to quickly identify promising drug candidates. The structural information of recently released crystal structures of main protease (Mpro) in APO and complex with inhibitors, N3, and 13b molecules was utilized to explore the binding site architecture through Molecular dynamics (MD) simulations. The stable state of Mpro was used to conduct extensive virtual screening of the aforementioned drug pool. Considering the recent success of HIV protease molecules, we also used anti-protease molecules for drug repurposing purposes. The identified top hits were further evaluated through MD simulations followed by the binding free energy calculations using MM-GBSA. Interestingly, in our screening, several promising drugs stand out as potential inhibitors of Mpro. However, based on control (N3 and 13b), we have identified six potential molecules, Leupeptin Hemisulphate, Pepstatin A, Nelfinavir, Birinapant, Lypression and Octreotide which have shown the reasonably significant MM-GBSA score. Further insight shows that the molecules form stable interactions with hot-spot residues, that are mainly conserved and can be targeted for structure- and pharmacophore-based designing. The pharmacokinetic annotations and therapeutic importance have suggested that these molecules possess drug-like properties and pave their way for in-vitro studies.Communicated by Ramaswamy H. Sarma.

A Review on the Progress and Prospects of Dengue Drug Discovery Targeting NS5 RNA- Dependent RNA Polymerase.

Dengue virus (DENV) infection threatens the health and wellbeing of almost 100 million people in the world. Vectored by mosquitoes, DENV may cause a severe disease in human hosts called Dengue hemorrhagic fever (DHF)/Dengue shock syndrome (DSS), which is not preventable by any known drug. In the absence of a universally-accepted vaccine, a drug capable of inhibiting DENV multiplication is an urgent and unmet clinical need. Here we summarize inhibitory strategies by targeting either host biochemical pathways or virus-encoded proteins. A variety of approaches have been generated to design Directly-acting anti-virals or DAAs targeting different DENV proteins, with diverse success. Among them, DAAs targeting genome replicating viral enzymes have proven effective against many viruses including, Human Immuno-deficiency Virus and Hepatitis C Virus. DAAs may be derived either from existing compound libraries of novel molecules and plant secondary metabolites or devised through Computer-aided Drug design (CADD) methods. Here, we focus on compounds with reported DAA-activity against the DENV RNA-dependent RNA polymerase (RdRp), which replicate the viral RNA genome. The structure-activity relationship (SAR) and toxicity of the natural compounds, including secondary plant metabolites, have been discussed in detail. We have also tabulated novel compounds with known anti-RdRp activity. We concluded with a list of DAAs for which a co-crystal structure with RdRp is reported. Promising hit compounds are often discarded due to poor selectivity or unsuitable pharmacokinetics. We hope this review will provide a useful reference for further studies on the development of an anti-DENV drug.

Insights into structural dynamics of allosteric binding sites in HCV RNA-dependent RNA polymerase.

Inhibition of the viral RNA-dependent RNA polymerase (RdRp) to resolve chronic infection is a useful therapeutic strategy against Hepatitis C virus (HCV). Non-nucleoside inhibitors (NNIs) of RdRp are small molecules that bind tightly with allosteric sites on the enzyme, thereby inhibiting polymerase activity. A large number of crystal structures (176) were studied to establish the structure-activity relationship along with the mechanism of inhibition and resistance between HCV RdRp and NNIs at different allosteric sites. The structure and the associated dynamics are the blueprint to understand the function of the protein. We have implemented the ligand-based pharmacophore and molecular dynamic simulations to extract the possible local and global characteristics of RdRp upon NNI binding and the structural-dynamical features possessed by the known actives. Our results suggest that the NNI binding induces significant fluctuations at the atomic level which are critical for enzymatic activity, with minimal global structural alterations. Residue-wise mapping of interactions of NNIs at different sites exhibited some conserved interaction patterns of key amino acids and water molecules. Here, the structural insights are explored to understand the correlation between the dynamics of protein's subdomains and function at the molecular level, useful for genotype-specific rational designing of NNIs.Communicated by Ramaswamy H. Sarma.

Conformational Characterization of Linker Revealed the Mechanism of Cavity Formation by 227G in BVDV RDRP.

RNA-dependent RNA polymerase (RdRp) is a relevant antiviral drug target. We investigated a potent benzimidazole inhibitor (227G; IC50 = 0.002 µM) against bovine viral diarrhea virus (BVDV) RdRp; however, its inhibition action was completely impaired in the presence of a resistant mutation, I261M. The binding of 227G in mutant RdRp affected the binding site loop conformations (especially Linker) that increased the volume of the binding site. It was also observed that the innate Linker's flexibility was retained, which was otherwise completely frozen in the wild-type complex. The functional role of Linker was hypothesized that it is a multidocking site for RNA template, inhibitors, and the other proteins involved in replication complex formation. The binding phenomenon requires significant molecular flexibility and the large-amplitude conformational dynamics of Linker, which is currently unknown. We observed a bidirectional "hinge"-like motion of Linker from crystal position, indicating its pronounced flexible behavior. This study underscores the importance of Linker's flexibility in the functionality of BVDV RdRp and proposes the template entrance site for selective anti-BVDV drug discovery.